Possible effect of collective modes in zero magnetic field transport in an electron-hole 

bilayer 
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We report single layer resistivities of 2-dimensional electron and hole gases in an electron-hole 
bilayer with a lOnm barrier. In a regime where the interlayer interaction is stronger than the 
intralayer interaction, we find that an insulating state (dp/dT < 0) emerges at T ~ 1.5K or lower, 
when both the layers are simultaneously present. This happens deep in the "metallic" regime, even 
in layers with fc_p£ > 500, thus making conventional mechanisms of localisation due to disorder 
improbable. We suggest that this insulating state may be due to a charge density wave phase, 
as has been expected in electron-hole bilayers from the Singwi-Tosi-Land-Sjolander approximation 
based calculations of L. Liu et al [Phys. Rev. B, 53, 7923 (1996)]. Our results are also in qualitative 
agreement with recent Path-Integral-Monte-Carlo simulations of a two component plasma in the 
low temperature regime [ P. Ludwig et al. Contrib. Plasma Physics 47, No. 4-5, 335 (2007)] 
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INTRODUCTION 



The total energy of a system of electrons can be 
thought of as the sum total of the kinetic energies of the 
free particles and their potential energies due to mutual 
Coulomb interaction. The relative importance of the two 
contributions is measured by the parameter r s — E ee /Ef 
(where E ee = e 2 \j '(it N) / '4tt ne and Ef — nh 2 N/m e ff 
in 2-dimensions, with N electrons per unit area). r s is a 
crucial parameter (apart from disorder) that governs the 
behaviour of the system to a large extent. For example 
it is known that a system of electrons or holes would 
be "gas like" at r s w 1, "liquid like" at r s ss 10 and 
possibly a solid "Wigner crystal" at r s w 100, provided 
disorder driven localisation does not dominate and 
pre-empt the observable effects of interaction. Confining 
a large number of particles in a small area makes the 
inter-particle spacing small and hence the Coulomb 
repulsion large, but the kinetic energy of the particles 
increases even faster - making r s smaller. This somewhat 
counter-intuitive fact is a straightforward consequence 
of Fermi statistics and is true in all dimensions. 
Let us now consider two parallel layers of electrons or 
holes with 10 11 cm -2 electrons in each. As they are 
brought closer to each other the particles in one layer 
not only interact with others in the same layer but 
also with those in the other layer. The inter-particle 
spacing in the same layer stays fixed and is about 30nm. 
It is possible (though highly non-trivial) to make the 
distance between the two layers about lOnm without 
their wavefunctions beginning to overlap. lOnm is 
approximately the excitonic Bohr radius in Gallium 
Arsenide (GaAs) and is an important length scale. We 
thus get an electron to "see" another electron (or hole) 
only lOnm away, without paying the kinetic (Fermi) 
energy cost, because the two layers continue to be 
two separate Fermi systems. As a consequence bilayer 
systems can give rise to interaction-driven phases that 



are not possible in single layers. The case of attractive 
interaction (electron-hole) is more interesting. A rich 
phase diagram is anticipated for the ground state of a 
spatially separated electron-hole bilayer (EHBL) that in- 
cludes a superfluid 1 "^^, charge density waves (CDW), 
Wigner crystals (WC)2»Z£, an excitonic supersolid^ and 
a possible crossover from a Bose-condensate to Bardeen- 
Cooper-Schreiffer type stated. Recent techniques for 
making independent ohmic contacts to 2-dimensional 
electron gases (2DEG) and 2-dimensional hole gases 
(2DHG) 10-20nm apart have enabled transport measure- 
ments down to millikelvin temperature o 11 ' 12 ' 13 . This has 
greatly extended the range of densities and temperatures 
that could be explored in the first attempt to measure 
transport in EHBLs 1 ^. 

Quite naturally the initial experiments were focussed 
on measuring the Coulomb drag in these bilayer o 15 ' 16 . 
These results strongly suggest the emergence of a non 
Fermi-liquid phase at temperatures below T ~ IK. 
The scattering rate between the electrons and holes 
(that is measured in Coulomb drag experiments) is 
not determined by the bare Coulomb potential, as 
the inter-particle interaction is screened and it is the 
dielectric function e(q,ui), in the presence of interactions 
(and possible pairing fluctuations 4 ) that determines 
the measured Coulomb drag. On the other hand the 
dominant contribution to the single layer resistance 
comes from the scattering due to residual charged 
background impurities/traps and the remote ionised 
dopants (if any). The bare Coulomb potential due to 
these is screened by the same dielectric function. Thus 
it is worth asking the question: Does the interlayer 
interaction also bring about some distinct feature in the 
single layer resistances? This is the central question we 
investigate experimentally in this paper. A connection 
between drag and single layer resistance was earlier 
investigated in the context of in-planc magnctoresistance 
and magnetodrag in hole gases 1 ^. There is however a 
difference between the screening of interlayer electron- 
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hole interactions and impurity potentials. The electrons 
and holes are mobile and the impurities are static. 
The potentials due to electrons/holes change with time 
and one needs to consider dynamic screening. On the 
other hand static screening e(q, u) — 0) is sufficient 
to account for the effect of ionised impurities. Thus 
the Coulomb drag is sensitive to more features in the 
screening function (such as plasmon modes). Also, 
it is easy to imagine situations where the Coulomb 
drag and single layer resistivities would have opposite 
temperature dependences. For example, if the layers 
are individually insulating, then the resistivity of the 
layer will increase as the temperature is decreased, but 
the Coulomb drag will go to zero monotonically at low 
temperatures. While we are aware of these distinctions 
between drag and single layer resistance, we find striking 
features in the temperature dependence of the single 
layer resistivities of both layers. These features are 
unambiguously attributable to the presence of the other 
layer and hence electron-hole interaction only. 



II. DEVICES AND EXPERIMENTAL 
METHODS 

In our experiments the electron and hole den- 
sities (n,p)were varied between 4 x 10 10 cm -2 and 
1.6 x 10 11 cm -2 , separated by a lOnm Alo.9Gao.1As 
barrier. In this regime the inter-layer interaction is 
comparable to or even stronger than the intra-layer 
interaction. Several important length scales in the 
problem are very close to each other. Considering 
a typical density of 1 x 10 11 cm -2 , the inverse Fermi 
wavevector is fc^ 1 = l/v / 2~7m «13nm. The hole quan- 
tum well is 20nm wide. The peak-to-peak distance 
of the wavefunctions is d ss 25nm. For the lowest 
density (n = p = 4 x 10 10 cm -2 ), we have kpd = 1.2. 
At the same time, the (conduction band) excitonic 
Bohr radius in GaAs is as ~10nm. The interaction 
parameter for the electrons, r s (ratio of the interparticle 
Coulomb energy to the Fermi energy) at the lowest 
density is 2.8 and is «12 for the holes. Thus, we 
reach a regime d ~ r s as where excitonic phases^and 
density modulated phases^ have been predicted. The 
Molecular Beam Epitaxy growth, fabrication process 
and operation of these devices have been described in 
detail earlie r 11 ' 13 , so we do not describe them here. Two 
samples (Device D & E) fabricated from the same wafer 
(ID:A4268) were used in this study and gave very simi- 
lar results. (A third sample also gave similar results.) 
Coulomb drag data from device D was reported earlier^. 
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FIG. 1: (Color online)Data from device E, measured in a 
pumped He 3 cryostat. (a) p xx (T) of the hole layer, when no 
electrons are present, (b) p xx (T) of the hole layer with elec- 
trons present. Exactly the same current and voltage contacts, 
measuring current (~ lOnA) and frequency (7Hz) were used 
for consistency, (c) Shubnikov de-Haas data corresponding 
to the traces in (a) showing normal 2DHG behaviour. The 
lowest trace corresponds to p — 9 x 10 10 cm 2 . (d) "Metal- 
lic" behaviour of the holes at p = 9 x 10 10 cm 2 , with n — 0. 
(e)Near exponential rise of p xx (T), at p — 9 x 10 10 cm 2 with 
n — p. Note that the figure is plotted in a log vs i scale. 



example 1 -, at the cost of sacrificing independent contacts 
and independent tuning of the layer densities. These ex- 
periments, if they can be done on EHBLs in future, would 
certainly provide valuable data. However we show that 
even dc transport in EHBLs (with its limitations), yields 
results strongly suggestive of the appearance of CDW or 
WC phases in EHBLs. 



The complexity of these samples, at present, stands in 
the way of integrating coplaner waveguides and perform- 
ing high frequency (~ 100 MHz) experiments to probe 
the "pinning modes" of density modulated phases. Such 
experiments have been done on bilayer hole gases, for 



III. RESULTS 

We first show the temperature dependence of the re- 
sistivity of the holes, when no electrons are present (Fig. 



3 



[TJl) . This was achieved by keeping the interlayer bias 
voltage below the threshold (V e h ~ 1.55V) for accumu- 
lation of electrons. The density of the holes was con- 
trolled by using the backgate. The data shows features 
typical of numerous studies on Si-MOSFETS and GaAs 
based single layer 2DEGs or 2DHGs in the context of 
the 2D Metal-Insulator transition (MIT) problem. At 
low densities (p<5 x 10 10 cm -2 ) the traces are clearly in- 
sulating (dp/dT < 0), and turn metallic (dp/dT > 0) 
at higher densities. The crossover happens as expected 
near p xx ss h/e 2 or cquivalcntly when kpl ss 1. The 
central result of this paper is the strikingly different be- 
haviour of the 2DHG when it sees the 2DEG. Fig. [Tp, 
where we plot the hole resistivity for several equal elec- 
tron and hole densities, shows that even the traces which 
were metallic in Fig. []Jl (e.g. p = 9 x 10 10 cm~ 2 in 
Fig- QJ) turn very clearly insulating by T ~ IK even at 
kpl > 50. We will show later that n — p is not essential 
for this insulating behaviour to appear. Simple argu- 
ments based on localisation due to background impuri- 
ties/dopants/defects cannot account for this because the 
2DHG (without the 2DEG) would see exactly the same 
charged background impurities. This is the key point 
that allows us to attribute this effect to electron-hole in- 
teractions. In fact if we look at the high temperature 
end of the traces (T ps 1.5K) we find that the resistivity 
of the 2DHG actually decreases significantly as the elec- 
trons are introduced. For the p = 6 x 10 10 cm~ 2 trace, 
p xx (n = 0) = 3808fi/n, whereas p xx (n = p) = 12030/n 
at T = 1.5K (compare Fig. QJa) and DJb)). This is 
a very significant (threefold) increase in mobility and 
consistent with the findings of Morath et al in EHBLs 
with 30nm barriers^. The effect is much stronger here, 
most likely because of the much thinner barrier used in 
this study. The increase in mobility at the high tem- 
perature end is most likely brought about by the extra 
screening of charged impurities due to the 2DEG and the 
"squeezing" of the hole wavefunction itself. The Coulomb 
drag was also measured for both the samples (data not 
shown) and at T = 1.5K the values obtained are consis- 
tent with a T 2 behaviour. Below T ps IK both samples 
show a strong deviation from T 2 behaviour, as reported 
recently by us and Seamons et a ^ 15 ' 16 . (The T 2 behaviour 
can be well understood within Fermi-liquid theories, but 
the strong non-monotonic deviations cannot be.) But a 
simple Mathiessen's rule based addition of two contri- 
butions due to Coulomb drag and impurity scattering 
cannot explain an increase in mobility. At the same time 
(for n=p=9xl0 10 cm -2 ) the measured Coulomb drag at 
T < IK is much smaller than the change in single layer re- 
sistivity (of the 2DHG) brought about by the presence of 
the 2DEG. Simple (incoherent) addition of various scat- 
tering rates no longer appears to be able to account for 
the observed magnitude of the effect. 

If we now consider the behaviour of the electron layer, 
at n — p = 9 x 10 10 cm -2 , the situation appears even 
more striking because the insulating behaviour can be 
easily seen at kpl > 500. Note that (in Fig. ^ although 
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FIG. 2: (Color online) A comparison of the electron(blue) 
and hole (red) resistivity data from device D, measured in a 
dilution refrigerator with base temperature of w 50mK. At 
this density were p e = 1.5 x 10 6 cm 2 V~ 1 s~ 1 & p,h = 1.1 x 
10 5 cm 2 V~ 1 s~ 1 (at T » 1.5K). The inset shows the amount 
of conductivity lost by the layers as a function of temperature. 

the electron and hole densities are same, the single layer 
resistances differ by nearly a factor of 10 at T = 1.5K 
and by a factor of 80 at T sa 50mK. However the amount 
of conductivity (Aa xx ) lost by the layers as the tem- 
perature is reduced below 2K, does not differ by more 
than a factor of 2. Aa xx between 2K and rs 50mK, is 
of the order of 2mlS, which places it much beyond what 
weak localisation based effects (Ac ~ e 2 /h 40/ilS) can 
account for—. In Fig. [3] we show how dp/dT of both 
the electron and hole layers below ~ 0.7K turn from 
"metallic" to "insulating" as the carrier density of the 
electron layer is increased in steps, while keeping the hole 
density (p — 1.6 x 10 cm" 2 ) constant. We chose this 
density, such that the sheet resistivities (at T ~ 1.5K) 
could be made nearly same for both the layers. The in- 
sulating state gradually evolves from a "metallic" state 
as the holes see more electrons. As far as we can say, the 
insulating state appears simultaneously in both layers. 
However, n = p does not appear to play a particularly 
special role in the problem, as has been foun d 15 ' 16 in the 
context of Coulomb drag in these bilayers. 



IV. DISCUSSION 

To understand what might be causing this we organise 
our discussion in three steps. First we have already 
mentioned several points which show that the physics of 
single layer 2DHG or 2DEG cannot explain the results. 
Next we consider whether any inhomogeneity caused 
by the large voltage bias across the lOnm barrier can 
affect our results. Finally we discuss the possibility of 
the emergence of collective modes (like CDWs) in an 
EHBL and refer to theoretical works that suggest that 
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FIG. 3: Data from device E, measured in a pumped He 3 cryo- 
stat. The hole density was kept fixed at p — 1.6 x 10 11 cm -2 
while the electron density was varied in steps n—4, 6, 8 x 
10 10 cm -2 . (a) The p xx (T) for holes is "metallic" when the 
number of electrons is n — 4 x 10 10 cm -2 on the other side of 
the barrier. As the electron density is increased, the p xx (T) of 
the holes becomes insulating, (b), (c) & (d) show the p xx (T) 
of the electrons. Note that the electrons are metallic at low 
density and insulating at the two higher densities, IcfI ~ 170 
or higher for all the traces. 



such modes are more likely in EHBLs than 2 x 2DEGs or 
2x2DHGs. 

Let us now consider the possible effect of some variation 
in thickness of the barrier, resulting in spatial variation 
of the interlayer capacitance. The question here is: 
In our devices does this lead to coexisting low density 
(insulating) regions and high density (metallic) regions? 
If we assume a variation of a monolayer («0.5nm), 
this would amount to a maximum of 5% fluctuation. 
The corresponding density fluctuations would not (for 
example) be able to force regions with p = 9 x 10 10 cm -2 
and regions with p«5 x 10 10 cm~ 2 (which is where we 
find the 2D MIT in the 2DHG when no 2DEG is present) 
in the same sample to coexist. Unless the inhomogeneity 
is so strong that the low density regions become strongly 
insulating, it cannot be the driving factor. If the 2DHG 
did become strongly inhomogeneous, it would have 
increased the sheet resistance at high temperatures as 
well. As long as the underlying idea is based on the 
general (single layer) picture of insulating behaviour 
above h/e 2 and metallic behaviour below this resistance, 
it is not possible to explain the results. A fundamentally 
different mechanism that can generate an insulating 
state at p xx « h/e 2 must be sought here. 
We look for an approach that specifically takes into ac- 
count that in a bilayer there can be collective modes that 
have no analogue in single layer situation. Collective 
modes are zeros of the dielectric function e(q,u>) where 
a system can develop density modulations without 
an external perturbation. For a bilayer e(q, u>) is a 
2x2 matrix, whose determinant would be zer o 21 ' 22 . 
Equivalently, the charge susceptibility of the system 



would diverge at these points. If such a mode occurs 
at u=0, it can be a CDW or WC phase, depending on 
the wavevector at which the divergence occurs. In a 
bilayer an electron "sees" another electron (or hole) only 
lOnm away, without paying a kinetic (Fermi) energy 
cost. As a consequence bilayer systems can give rise 
to interaction-driven phases more easily than single 
layers. In fact, existing calculations^ of the bilayers for 
particle densities corresponding to r s ~ 1 — 10 mention 
the possibility that a divergence in the eigenvalues of the 
bilayer e(q, u = 0) matrix is easier to get in an EHBL 
compared to 2x2DEGs or 2x2DHGs. This point needs 
attention because experimental data from 2x2DEGs 
or 2x2DHGs with a barrier separation of lOnm and 
densities as low as lxl0 lo cm -2 have been available 
for some tim e 23 ' 24 . These samples had similar electron 
mobilities as the EHBLs reported here, but higher hole 
mobilities. The principal interest with these devices has 



been focussed on the 
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bilayer Quantum Hall 



state, but interestingly no collective mode at B=0 has 
been reported in these^. Following Liu et al^ we con- 
sider the expression for eigenvalues of the susceptibility 

matrix, one of which must diverge to support a CDW. 
( \ _ 2 

~ Xe- 1 +X^ 1 ±V / (Xc- 1 -X^ 1 ) 2 +4{(l-Ge h ( 9 ))V eft ( 9 )} 2 ' 

Xe,h are the single layer responses to the external 

2ne 2 



exp(-qd) is the Fourier 



potential and V e h{q) 

transform of the interlayer Coulomb potential. The 
crucial role is played by the local field correction G e h(q), 
often calculated using the Singwi-Tosi-Land-Sjolandcr 
(STLS) approximation 26 that relates the local field to the 
structure factor. The interlayer local field correction is 
larger in magnitude in case of an attractive potential 6,8 . 
Due to this a divergence of the in-phase mode in an 
EHBL is easier to get than the corresponding divergence 
in the out-of-phase mode in a 2x2DEG or a 2x2DHG. 
At the same time the strong peak in x(o)j (that signals 
the formation of a CDW state) has been found to be 
strongly temperature dependent, in numerical studies^. 
The range of temperatures at which the rapid increase in 
resistance is seen is consistent with the expectation that 
it should be at a temperature smaller than the Fermi 
temperature of the layer with heavier mass (Tp w 5K 
for holes, at p = 9 x 10 10 cm~ 2 ). 

Path Integral Monte Carlo simulation of a spatially 
separated 2-component plasma has been carried out by 
Ludwig et for parameters relevant to GaAs-AlGaAs 
at low temperatures. The real space distribution of 
the holes (in presence of electrons) started developing 
density modulations, while the lighter electron liquid 
remained in a much more homogeneous state. Such a 
situation is quite consistent with our observations of an 
insulating behaviour of the holes, and a less pronounced 
(but clearly discernible) effect on the electrons. A 
density modulation (at the same wavevector) with 
larger amplitude in the hole layer and smaller amplitude 
in the electron layer is also supported by the STLS 
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calculations^. 

There is nothing in our observations that rules out a 
much awaited excitonic state in electron-hole bilayers. 
The interlayer pair correlation in an excitonic (bound) 
state would involve a more drastic modification of the 
structure factor and the local field corrections, compared 
to that of a CDW state. A CDW may well be a precursor 
to an excitonic stated 

Finally, we observe that a CDW may also be a possible 
mechanism behind the temperature dependent enhance- 
ment in Coulomb drag reported recentl y"-* - 6 . We con- 
sider a situation where the densities (n,p) in each layer 
undergo an in-phase modulation S with a long wavelength 
(i.e. small q/kp), such that the densities in successive 
half periods are n — S and n + S. We assume that we 
can evaluate the contribution to drag over each half cy- 
cle and add the voltages in series. Over each half cy- 
cle the contribution is still assumed to be proportional 



to l/(np) 3 / 2 , where n,p are now the local densities- 7 . 
Thus the contribution of two adjacent half cycles would 
be oc {n+s)3/2 \ p+s)3/2 + {n _ s)3/ 2 (p _ s)3/ 2 ■ Keeping terms 
upto 8 2 we get the contribution from each CDW period 

fn hp nc 2 l + (9/4),5 2 /np+(3/8)<5 2 (l/n 2 + l/p 2 ) ™ . 

De ^ (« P )3/ 2 (i-5 2 /n 2 )3/2(i_52 /p 2 ) 3 / 2 ■ i *iis quan- 
tity is clearly larger than the unperturbed contribution 
( np )3/2 1 an d it is clearly temperature dependent, since 
the amplitude of a CDW (i.e 5) can increase as the tem- 
perature decreases. This argument neither depends on 
the exact form of the power law nor the matching of the 
densities. The fact that it is of order 5 2 may have some 
significance in explaining why the observed effec t 15 ^ 16 is 
rather small and not a large peak or divergence expected 
at matched densities from an excitonic phase. 
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